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A chiral rhenium complex with predicted high parity violation effects: 
synthesis, stereochemical characterization by VCD spectroscopy and 
quantum chemical calculations 


Nidal Saleh, a Samia Zrig, a,& Thierry Roisnel, a Laure Guy, c Radovan Bast/ Trond Saue/ 
Benoit Darquie,^ and Jeanne Crassous* a 


With their rich electronic, vibrational, rotational and hyperfine structure, molecular systems have the potential to play a deci¬ 
sive role in precision tests of fundamental physics. For example, electroweak nuclear interactions should cause small energy 
differences between the two enantiomers of chiral molecules, a signature of parity symmetry breaking. Enantioenriched ox- 
orhenium(VII) complexes and R -(+)~3 bearing a chiral 2-methyl-1-thio-propanol ligand have been prepared as potential 
candidates for probing molecular parity violation effects via high resolution laser spectroscopy of the Re=0 stretching. Although 
the rhenium atom is not a stereogenic centre in itself, experimental vibrational circular dichroism (VCD) spectra revealed a 
surrounding chiral environment, evidenced by the Re=0 bond stretching mode signal. The calculated VCD spectrum of the R 
enantiomer confirmed the position of the sulfur atom cis to the methyl, as observed in the solid-state X-ray crystallographic struc¬ 
ture, and showed the presence of two conformers of comparable stability. Relativistic quantum chemistry calculations indicate 
that the vibrational shift between enantiomers due to parity violation is above the target sensitivity of an ultra-high resolution 


infrared spectroscopy experiment under active preparation. 

1 Introduction 

Chirality is a challenge in diverse domains of chemistry, such 
as in the development and synthesis of molecules for pharma¬ 
ceuticals, 1 agrochemicals, flavors, in nanotechnology 2 3 or in 
asymmetric catalysis. 4 Furthermore, studying unconventional 
chiral systems may be of great help for understanding the ori¬ 
gin of prebiotic homochirality. 5 More recently, very simple 
chiral molecules have been the object of interest for parity vi¬ 
olation (PV) measurements, aiming at probing a tiny energy 
difference between the right- and left-handed enantiomers of 
a chiral molecule. ^ 8 During the last decade, we have been in¬ 
terested in observing this effect using an ultra-high resolution 
infra-red (IR) spectroscopy technique on a supersonic molecu¬ 
lar beam. 913 This collaboration between theoretical chemists, 
synthetic chemists, and spectroscopists, aiming at observing a 
vibrational PV frequency shift between enantiomers of chiral 
molecules, has recently been reviewed. 12 
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A number of chiral compounds with heavy metals as the 
chiral centers have been theoretically studied for parity vi¬ 
olation. 14 19 In this context “chiral at metal” oxorhenium 
complexes have proven particularly suitable heavy transition 
metal complexes because they may fulfill all the requirements 
needed for a successful PV experiment. Indeed, the best can¬ 
didate for our experimental test should i) show a large PV vi¬ 
brational frequency difference of an intense fundamental band 
preferably within the CO 2 laser operating range (900-1100 
cm -1 ), ii) be available in large enantiomeric excess or, ide¬ 
ally, in enantiopure form and at the gram-scale, iii) have a 
reasonably simple structure so as to facilitate the spectroscopy 
and maintain a favourable partition function, and iv) allow the 
production of the supersonic expansion, thus sublime with¬ 
out decomposition. 12 Only a few simple “chiral at metal” rhe¬ 
nium complexes are known in the literature. 20 29 Two classes 
of chiral complexes have been studied for PV experiments 
(see Figure [I]): i) oxorhenium complexes such as 1, bearing 
hydrotris(l-pyrazolyl)borate (Tp) ligand and a chiral bidentate 
ligand, 23 and ii) enantiopure “3+1” oxorhenium complexes 
bearing a tridentate sulfurated ligand and a monodentate halo¬ 
gen (such as 2) or a chalcogenated ligand. 24 25 Methyltriox- 
orhenium (MTO) is an efficient catalyst for reactions such as 
for instance epoxidation of olefins J 30 More interestingly for 
our purposes, it sublimes very easily 13 (vapour pressure of a 
few tenths of mbars at room temperature). Furthermore, a few 
examples of sublimable MTO derivatives have been described 
in the literature 29 and the strategy of preparing simple chiral 
easily sublimable MTO derivatives seems an appealing route 





to find a suitable candidate molecule for the experimental PV 
test. Our target molecule was therefore the enantiopure ox- 
orhenium complex 3 which is a chiral derivative of MTO bear¬ 
ing a 2-methyl-thiopropanolate ligand. In this article we detail 
the synthesis, the stereochemical characterization and vibra¬ 
tional circular dichroism (VCD) spectroscopy of complex 3 
to examine the dissymmetric environment around the rhenium 
atom. We furthermore carry out relativistic quantum chemical 
calculations of the vibrational frequency shift associated with 
the parity violation energy difference between enantiomers of 
complex 3. Experimental and computational details are given 
in sections [ 7 ] and [ 8 ] respectively. 
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Scheme 1. Synthesis of the enantioenriched ligand 
2-methyl-1-thiopropanol R-(-)-4 from R -(+)-propylene oxide. 
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Fig. 1 Examples of chiral oxorhenium complexes studied for parity 
violation effects (refs. [23ll25l and this work). 


2 Synthesis 

The synthesis of the enantiomers of the 2-methyl-thiopropanol 
(R-(-)- and S-(+)~ 4) ligands was achieved using enantiomer- 
ically pure propylene oxide as the starting precursor as de¬ 
scribed in Scheme 1. The regioselective opening of R-(+)~ and 
^-(-^propylene oxide with trityl thiol 31 gave respectively the 
alcohols R-(+)~ and 5-(-)-5 with 90% yield. The subsequent 
deprotection followed by in situ oxidation with iodine yielded 
quantitatively disulfide R,R-(-)~ and S,S-(+)-6, which were 
finally reduced by LiAlH 4 to R-(-)- and S-(+)-4 with 84% 
yield. Their optical rotation and all spectroscopic character¬ 
istics are in agreement with already published results.^ 32 The 
(+)- and (-)-3 enantiomeric complexes were then prepared 
by reacting the enantiopure R-(+)~ and 4 with MTO in 
dry dichloromethane according to a previously reported pro¬ 
cedure 33 as depicted in Scheme 2. 

3 X-ray crystallography 

Single crystals of a racemic sample of the oxorhenium(VII) 
complex 3 were obtained after sublimation under vacuum. 
Due to significant disorder of the methyl group coming from 
the two possible configurations of the stereogenic carbon cen¬ 
ter, the X-ray structure of complex 3 was solved in the non- 
centrosymmetric P21 space-group. A final Flack parameter of 
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Scheme 2 Synthesis of the enantioenriched oxorhenium complexes 
R-(+)- and from MTO and enantioenriched 

2-methyl-1-thiopropanol R-(-)~ and S-(+)-4. 


0.46 was obtained, indicating the presence of a racemic mix¬ 
ture. As depicted in Figure [2] the rhenium center is penta- 
coordinated, the Cl atom is trans to the 03 oxygen, and the 
Re,03,S,Cl lie in the same plane, while the two oxo groups 
are symmetrically placed on each side of this plane (for exam¬ 
ple, angles 03-Re-01: 100.66° and 03-Re-02: 105.03° for 
one of the four independent molecules of the unit cell). Con¬ 
sequently the oxygen atoms 01 and 02 only differ in their 
chemical environment due to the proximity of the asymmetric 
carbon C2. The Re-0 bond lengths in the two oxo groups are 
classic (1.675 A to 1.687 A). Since oxorhenium complexes 
most often crystallize in a square pyramidal geometry, with 
an oxygen atom of the Re=0 bond being placed in the apical 
position, the geometry around the rhenium atom is rather un¬ 
common. 34 Finally, in the solid state the methyl group C3 is 
placed in the equatorial position, and the dihedral angle 03- 
C2-C4-S within the five-membered ring is about -51°. The 
conformational space of complex 3 has been further explored 









by density functional theory (DFT) calculations, as discussed 
in the next section. 



Fig. 2 X-ray crystallographic structure of the complex R-(+)~ 3. 


4 VCD spectroscopy 

For our fundamental physics purposes we would like simple 
chiral molecules containing one or more heavy metals, prefer¬ 
ably as a stereogenic center or at least in its vicinity. The for¬ 
mer is nonetheless not straightforwardly realized.*^ 25 From 
the experimental viewpoint, VCD spectroscopy, which has 
been under-explored in chiral complexes, can be a very pow¬ 
erful tool to study the dissymmetric environment around the 
metal center. 35 38 In order to get more insight into the latter, 
the IR and VCD spectra of complex 3 were calculated and 
compared to experiment. 

Our first computational goal, however, was to find a set of 
stable conformers of complex 3 for a subsequent calculation 
of IR and VCD spectra. Using the R enantiomer and starting 
from the X-ray structure in which the S-Re bond is placed 
in cis position to the Re-Me one, we carried out a potential 
energy scan along the 03-C2-C4-S dihedral angle at the DFT 
level using the B3LYP functional (def2-TZVPP basis^S). We 
identified two stable conformers, 3a-cl and 3a-c2 at the 03- 
C2-C4-S dihedral angles -39.8° and +37.5°, respectively. As 
shown in Figure [3j the conformers are separated by a barrier 
of almost 10 kJ/mol and differ in energy by 3.4 kJ/mol, the 
latter corresponding to Boltzmann populations of 80%:20% at 
296 K. 

The conformers 3a-cl and 3a-c2 are not the only possible 
realizations of the ligation and we have extended the study by 
optimizing the structures of corresponding 3b-cl and 3b-c2 
conformers where the O-CHMe-CF^-S ligand is connected 
to Re with S-Re trans to Re-Me. The conformers 3b-cl and 
3b-c2 are separated by 5.4 kJ/mol, corresponding to Boltz¬ 
mann populations of 90%: 10% at 296 K. Conformer 3b-cl is, 
however, 18.4 kJ/mol higher in energy compared to 3a-cl, and 
we will in the following be able to confirm the predominance 
of form 3a-cl from simulation of the spectra. Before turning 
to the discussion of simulated spectra we would first like to 
discuss the energetic separation between the four studied con¬ 
formers. To clarify whether the energetic preference for 3a-cl 
was sterically controlled by the Me group at the C2 atom or 



Fig. 3 Potential energy scan of complex 3 along the 03-C2-C4-S 
dihedral angle (B3LYP, def2-TZVPP basis, all other coordinates 
relaxed). 
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Fig. 4 B3LYP (def2-TZVPP basis) equilibrium structures of the 
theoretically studied conformers. 







electronically controlled according to the HSAB principle^ 
we have optimized structures for 3a’-cl, 3a’-c2, 3b’-cl, and 
3b’-c2 conformers, where the Me group at the C2 atom has 
been replaced by a H atom. By doing this we find that the con¬ 
formations cl and c2 are energetically almost equal for both 
3a’ and 3b’ (separated by less than 0.02 kJ/mol, but retaining 
the barrier heights) while the difference between 3a’-cl and 
3b’-cl remains 18.0 kJ/mol, as compared to the difference of 
18.4 kJ/mol between 3a-cl and 3b-cl. This indicates that the 
energy difference between cl and c2 is sterically controlled 
by the Me group attached to the C2 atom and can be tuned by 
bulkier ligands, while the preference for 3a over 3b is elec¬ 
tronically controlled by the O (HSAB harder) and S (HSAB 
softer) atoms.This theoretical study nicely illustrates the 
concept of the trans-effect in coordination chemistry. 

In order to verify the findings obtained using the relative en¬ 
ergies and to obtain more insight into the chiral environment 
around the metal center, we next compared simulated IR and 
VCD spectra of complexes 3a and 3b to the experimentally 
recorded spectra. The experimentally recorded IR spectrum of 
complex 3a in CD 2 CI 2 is shown as a dashed line in Figure [5] 
The IR signature is identical for both enantiomers (at this res¬ 
olution). The experimental VCD spectra of R -(+)~3 and S-(- 
)-3 enantiomers in CD 2 CI 2 are shown in Figure [6] in red and 
blue color, respectively. The experimental VCD spectra show 
a mirror-image relationship and exhibit a strong absorption 
band at 1012 cm -1 (dissymmetry factor Ae/e = 1.5 x 10 -4 ) 
which corresponds to the oxorhenium symmetric bond stretch¬ 
ing (the antisymmetric stretch, although not accessible to our 
VCD spectrometer, is also visible on the experimental IR 
spectrum of Figure[5]at ~ 980 cm -1 ). The presence of such a 
band indicates a chiral environment around the rhenium atom, 
although the rhenium is not a stereogenic center in itself. 

Comparing the observed IR and VCD spectra with the cal¬ 
culated spectra for the Boltzmann-averaged 3a and 3b we note 
that the match for conformer 3a is not perfect, but it is better 
than for 3b, which we expect based on energetic arguments. 
We should also note that the agreement between experimen¬ 
tal and simulated spectra seems to be better for the Stuttgart 
ECP/6-31Gd basis, already employed in Refs. |23]and|24l than 
for the def2-TZVPP basis, although the latter is a more flexi¬ 
ble basis which we consider close to basis set limit at the DFT 
level. This can probably be attributed to error cancellation ef¬ 
fects. Overall, the VCD spectrum calculated for the R- stere¬ 
ochemistry using the Stuttgart ECP/6-31Gd basis set repro¬ 
duces the experimental spectrum for the (+) enantiomer well, 
except for the peak at 1160 cm" 1 . The discrepancy may be due 
to the flexibility of the five-membered rhenacycle in solution 
that induces additional conformations. In conclusion this con¬ 
formational analysis and IR/VCD studies show that although 
the Re atom is not strictly speaking a stereogenic center, the 
Re-0 stretching bond at 1012 cm' 1 has a strong well-defined 



Fig. 5 Overlay of the simulated Boltzmann-averaged IR spectra for 
complex 3 (black solid lines, B3LYP, scaled), and the experimental 
spectrum (blue dashed line, 90 mg in 1 mL of CD 2 CI 2 , 100 ju m cell, 

4 cm -1 resolution, 3000 scans). 

chiroptical signature that is a testimony of the chiral environ¬ 
ment around the metal center. In our opinion this makes the 
enantiopure chiral oxorhenium(VII) complex 3 a good candi¬ 
date for a PV measurement. For this reason, the PV vibra¬ 
tional shifts associated with the Re-0 stretching modes were 
predicted by relativistic DFT calculations, as described in the 
following. 

5 Calculated PV shifts 


Table [T] summarizes the results of our computational study of 
the selected two conformers 3a-cl and 3a-c2 (see Figure [4]), 
where we have studied the symmetric and anti-symmetric 
Re=0 stretching mode frequencies using HF and two density 
functional approximations, the hybrid B3LYP and the non¬ 
hybrid PBE exchange-correlation functionals. Comparing the 
4-component DC Hamiltonian PV fundamental transition fre¬ 
quency differences Avo^i (Table [I]), we can first note that our 
calculations predict PV vibrational frequency differences in 
the sub-Hz range, above the anticipated experimental uncer¬ 
tainty in the frequency difference measurement. 12 This is an 
encouraging result for our collaboration. However, from the 





















Table 1 Calculated reduced masses (in amu), vibrational transition frequencies (in cm -1 ) of the anti-symmetric and symmetric Re=0 
stretches for the conformers 3a-cl and 3a-c2, with corresponding 4-component DC Hamiltonian PV fundamental transition frequency 
differences Avq^i (in Hz). 


Conformer 

Mode 

Method 

Red. mass 

Harmonic 

Fundamental 

PV shift 

3a-cl 

asym 

HF 

14.2074 

1106 

1102 

-0.211 



B3LYP 

15.7268 

986 

982 

0.078 



PBE 

15.8522 

951 

944 

0.195 


sym 

HF 

8.7326 

1174 

1165 

0.344 



B3LYP 

12.8095 

1022 

1010 

0.150 



PBE 

14.8015 

980 

968 

0.036 

3a-c2 

asym 

HF 

15.2424 

1106 

1102 

0.219 



B3LYP 

15.7515 

985 

981 

-0.119 



PBE 

15.8566 

950 

944 

-0.271 


sym 

HF 

8.4859 

1175 

1165 

-0.056 



B3LYP 

7.2491 

1021 

1010 

-0.170 



PBE 

13.1114 

979 

966 

-0.084 


theoretical point of view, the significant variation of the calcu¬ 
lated PV transition frequency differences is troubling, show¬ 
ing high sensitivity to both the method applied (HF or density 
functional approximations) and to structural changes. Com¬ 
paring HF with B3LYP and PBE we see that the order of mag¬ 
nitude can change as well as the sign. In principle we should 
rather trust the DFT calculations, since they incorporate elec¬ 
tron correlation, but the performance of DFT functionals with 
respect to PV shifts for these heavy element compounds needs 
further study (in progress in our laboratory). If we single out 
one method then a slight modification of the structural periph¬ 
ery seems to have a major effect on the chiral environment 
around the heavy atom (as exemplified by the sign change be¬ 
tween the PV shifts for 3a-cl and 3a-c2 in Table [TJ. This 
makes parity violation chemically interesting, as detailed anal¬ 
ysis of the underlying mechanism may provide a deeper un¬ 
derstanding of the electronic structure of chiral molecules. It 
should be emphasized that conformational averaging, as car¬ 
ried out in the simulation of IR and VCD spectra in section |4j 
is not an issue in the ultra-high resolution experiment since it 
will resolve individual lines of each conformed 

6 Conclusion 

We have reported the synthesis and resolution, characteriza¬ 
tion, and computational study of a new oxorhenium(VII) chi¬ 


ral compound. It is designed for the molecular beam exper¬ 
iment, proposed in Ref. Jj2i dedicated to the observation of 
PV effects via high-resolution vibrational spectroscopy. As an 
MTO derivative, this solid-state complex is expected to sub¬ 
lime easily, a key aspect to allow gas phase studies. Mea¬ 
surements and calculations show that the normal mode asso¬ 
ciated with the Re=0 stretch lies within the CO 2 laser oper¬ 
ating range. Enantioselective synthesis was developed for the 
preparation of enantioenriched complexes. DFT calculations 
allowed the identification of four conformers (3a-cl, 3a-c2, 
3b-cl and 3b-c2, see Figure]?]) as well as the rationalization of 
their relative stability. We find that the stability 3a with respect 
to 3b can be attributed to a trans effect, whereas the smaller 
energy difference of cl with respect to c2 arises from a steric 
effect. The experimental VCD spectra are accordingly dom¬ 
inated by the two most stable conformers 3a-cl and 3a-c2, 
as verified by their theoretical simulation which furthermore 
allowed determination of their absolute configuration. 

The VCD study has brought to light the chiral environment 
surrounding the rhenium atom, albeit not a stereogenic centre 
in itself. Motivated by this observation, the vibrational tran¬ 
sition frequency differences between the enantiomers due to 
PV was calculated using Hartree-Fock (HF) and DFT. These 
results confirm the 0.1 Hz order of magnitude for DFT PV vi¬ 
brational frequency differences, as reported earlier for oxorhe- 
nium complexes, which would be above the 0.01 Hz experi- 










Fig. 6 Overlay of the simulated Boltzmann-averaged VCD spectra 
for complex R -3 (black solid lines, B3LYP, scaled), and the 
experimental spectra (red and blue dashed lines correspond 
respectively to (+)-3 and (-)-3, 90 mg in 1 mL of CD 2 CI 2 , 100 jum 
cell, 4 cm -1 resolution, 3000 scans). 


mental uncertainty expected in the frequency difference mea¬ 
surement. ■ElThe present study puts emphasis on the sensitivity 
of the PV vibrational frequency difference to the chemical en¬ 
vironment around the rhenium centre, as reported earlier by 
us £ 25 This work represents an important step towards the first 
experimental observation of PV in molecular systems. 

7 Experimental details 

7.1 General 

Most experiments were performed using standard Schlenk 
techniques. Solvents were freshly distilled under argon from 
sodium/benzophenone (THF) or from phosphorus pentoxide 
(CH 2 CI 2 ). Starting materials were purchased from ACBR 
(MTO) or from Aldrich. The synthesis of ligand 4 enan¬ 
tiomers was performed by using a modified procedure . 32 Col¬ 
umn chromatography purifications were performed in air over 
silica gel (Merck Geduran 60, 0.063-0.200 mm). : H and 13 C 
NMR spectra were recorded on a Bruker AM300. 13 C NMR 
spectra at 50.4 MHz were recorded on a Bruker DPX 200. 
Chemical shifts were reported in parts per million (ppm) rel¬ 
ative to Si(CH 3)4 as external standard and compared to the 
literature. IR and VCD spectra were recorded on a Jasco 
FSV-6000 spectrometer. Specific rotations (in deg cm 2 g _1 ) 
were measured in a 10 cm thermostated quartz cell on a Jasco 
P1010 polarimeter. Elemental analyses were performed by 
the group CRMPO, University of Rennes 1. CCDC refer¬ 
ence number 910638 contains the crystallographic data for 
complex 3. These data can be obtained free of charge at 
www.ccdc.cam.ac.uk/conts/retreving.html or from the Cam¬ 
bridge Crystallographic Data Center, 12 union Road, Cam¬ 
bridge CB2 1EZ, UK; Fax: (internat.) +44-1223-336-033; E- 
mail: deposit@ccdc.cam.ac.uk. 

S-(-)-l-(tritylthio)propan-2-ol (S-(-)-5). ^r-BuLi (4.78 
mmol, 2.5 M, 1.91 mL) was added dropwise to a triphenyl- 
methanethiol solution (4.67 mmol, 1.29 g) in 15 mL dis¬ 
tilled THF cooled at 0°C, where a red-rose color persisted. 
S-(-)-propylene oxide (4.28 mmol, 0.3 mL) was then added 
dropwise at 0°C with a change in color to pale yellow. The 
reaction mixture was stirred over night, quenched with 20 % 
AcOH in 20 mL methanol, diluted with water, and then ex¬ 
tracted with ethyl acetate. Purification over silica gel column 
chromatography (pentane/ethyl acetate, 9:1) provided S-(-)-5 
with 90% yield. l U NMR (300 MHz, CDC1 3 ) 5 7.40-7.50 (6 
H, m), 7.19-7.37 (9 H, m), 3.43 (1 H, sxt, J = 6.2 Hz), 2.40 
(1 H, d, J = LI Hz), 2.38 (1 H, s), 1.07 (3 H, d, J = 6.1 Hz). 
[a ] 23 = -54.4 (C = 10 " 3 M, CH 2 C1 2 ) ref. [32| [a ] 23 = -27.5 (C 
= 3.35 x 10 " 3 M, CHCI 3 ). The same procedure was used for 
the preparation of the R-(+)-5 enantiomer [a ]^ 3 = +54.4 (C = 
10 " 3 M, CH 2 C1 2 ). 

5,S'-(+)-l,l’-disulfanediylbis(propan-2-ol) ((S,S)-(+)-6). 














S'-l-(tritylthio)propan-2-ol (S-(-)-5) (300 mg, 0.897 mmol) 
dissolved in 50 mL DCM/Methanol (9:1) solution was added 
in portions (30 min) over Iodine solution, 1.1 g in 500 mL 
DCM/Methanol (9:1). The reaction mixture was stirred at 
room temperature for 30 min, quenched with 10% aqueous 
Sodium thiosulfate and washed with brine. The organic layer 
was separated and the aqueous layer was extracted with ethyl 
acetate twice, dried over MgSCL, and concentrated in vacuum 
to provide a dark brown precipitate which was purified by sil¬ 
ica gel chromatography (5% ethanol/chloroform) to provide 
(25,2’5)-(+)-6 as a yellow-brown oil (99%). 'H NMR (300 
MHz, CDC1 3 ) 5 3.93-125 (1 H, m), 2.90 (1 H, dd, J = 13.6, 
3.4 Hz), 2.69 (1 H, dd, J = 13.7, 8.5 Hz), 2.41 (1 H, br. s.), 1.31 
(3 H, d, J = 6.0 Hz). 13 C NMR (75 MHz, CDCI 3 ) 8 65.9,47.5, 
22.0. [a]^ = +232 (C = 2.7 x 10 “ 3 M, CH 2 C1 2 ), ref. M [a]« 
= +199.3 (C = 1.96 x 10 ‘ 3 M, CHCI 3 ). The same procedure 
was used for the preparation of the enantiomer 

[ajg 5 = -236 (C = 2.7 x 10 ' 3 M, CH 2 C1 2 ). 

S-(+)-2-methyl-thiopropanol (£-(+)-4). 100 mg of ( S,S)~ 
l,l’-disulfanediylbis(propan-2-ol) 6 ) inTHF was added 
dropwise to a stirring solution of LAH (3 eq.) in THF at 0°C, 
then stirred for 24 hrs at 50° C. The reaction was quenched 
with AcOH and then extracted with ether, dried over MgSCL, 
and concentrated under vacuum to obtain £-(+)-4 (42.2 mg, 
84%) as a colorless oil. 'H NMR (300 MHz, CDCI 3 ) 8 3.77- 
3.86 (m, 1 H), 2.74 (dd, 1 H, J = 14.1, 3.6 Hz), 2.5 (dd, 1 H, J 
= 13.4 Hz, 7.4 Hz), 1.27 (d, 3 H, J = 6.0 Hz), [affi = +172 (C 
= 5.10 x 10 " 2 M, CH 2 CI 2 ). The same procedure was used for 
the preparation of the R-(-)-4 enantiomer, [a ]^ 3 = -166 (C = 
5.10 x 10 ' 2 M, CH 2 C1 2 ). 13 C NMR (75 MHz, CDC1 3 ) 5 68.3, 
32.8, 21.6P 

S-(-)-3 complex. Methyltrioxorhenium (110 mg, 0.44 
mmol) was dissolved in 10 mL distilled DCM, S-(+)-2- 
methyl-thiopropanol (0.44 mmol, 40.1 mg) was then added, 
and the reaction mixture was stirred overnight under argon. 
Solvent was removed under reduced pressure to obtain 
3 as a yellow-orange precipitate. : H NMR (300 MHz, CDCI 3 ) 
5 5.30 (1 H, sxt, J = 6.1 Hz), 3.84 (1 H, dd, J = 11.3, 5.3 Hz), 
3.46 (1 H, dd, J = 11.7, 7.2 Hz), 2.50 (3 H, s), 1.51 (3 H, d, J = 

6.1 Hz). 13 C NMR (75 MHz, CDC1 3 ) 8 90.6, 46.2, 29.7, 20.8. 
[a] 23 = -17 (C = 3.1 x 10 " 3 M, CH 2 C1 2 ). The same procedure 
was used for the preparation of the R-(+)-3 enantiomer and for 
racemic 3. [a] 23 = +17 (C = 3.1 x 10' 3 M, CH 2 C1 2 ). 

7.2 VCD measurements 

Samples of (+)-3 and (-)-3, « 9 mg /100 juh CD 2 CI 2 , were 
placed in a 100-/rm path length cell with BaF 2 windows. IR 
and VCD spectra were recorded on a Jasco FVS-6000 VCD 
spectrometer with 3000 scans acquired and averaged at 4 
cm ' 1 resolution. An overlay of the observed spectra for the 
two enantiomers, combined with the simulated Boltzmann- 


averaged VCD spectrum for complex R- 3 in two different ba¬ 
sis sets, is presented in Figure [ 6 ] The average VCD spectrum 
of the two enantiomers was used as the VCD baseline. 

8 Computational details 

The optimized molecular structures and the corresponding IR 
and VCD spectra were obtained using the Gaussian 09 pack¬ 
age 41 employing the default harmonic approximation 42 and 
using Kohn-Sham DFT together with the hybrid functional 
B3 lypE 1HI (employing the VWN303 correlation functional). 
In addition we have also employed HF and the PBE® func¬ 
tional to obtain the optimized structures and the harmonic 
force field for the subsequent Lpy calculations. We have 
used two different basis sets: Aldrichs’ def2- TZVPpEsMI (all 
atoms) to obtain a result close to the basis set limit, and the less 
flexible 6-3 lGd 4 ® (all atoms except Re) in combination with 
the Stuttgart/Dresden ECP-60-MWB 47 (Re) for compatibility 
with a previous VCD study of a similar compound . 24 The vi¬ 
brational frequencies were scaled by 0.97 and the calculated 
intensities were converted to Lorentzian bands with 4 cm -1 
half-width at half-maximum for comparison to experiment. 

In order to calculate the PV vibrational frequency shift we 
have performed single-point energy and PV energy calcula¬ 
tions using the DIRAC 12 program 49 ^ following the normal 
mode coordinates taken from the corresponding Gaussian 09 
calculation. The single-point energies, now based on the 
4-component relativistic Dirac-Coulomb (DC) Hamiltonian, 
calculated along the effective core potential normal mode pro¬ 
vided the anharmonic potential for the numerical solution of 
the vibrational wave functions by the Numerov-Cooley algo¬ 
rithm . 50 51 

The PV vibrational frequency shift v m ^ n for the m -+ n vi¬ 
brational transition (in our case 0 -+ 1 ) is obtained using 

hAv m ^ n = 2((n\E FW (q R )\n) - (m\E FW (q R )\m )), (1) 

where | n) and | m) are the corresponding Numerov-Cooley vi¬ 
brational wave functions. Epy(q R ) is the PV energy along the 
normal mode coordinate q of the R enantiomer, evaluated as 
an expectation value of the operator 

#PV = ^/7pv.x; (2) 

K 2v2 i 

in which appears the weak nuclear charge Q w ,k = Zk(1 — 
4sin 2 0w) — Nk with Z R and N R representing the number of 
protons and neutrons in nucleus K , respectively, and sin 2 0 w = 
0.2319 is the employed Weinberg parameter. The normalized 
nuclear charge densities p R , in this work represented by Gaus¬ 
sian distributions , 52 restrict integration over electron coordi¬ 
nates r i to nuclear regions, thus providing a natural partition¬ 
ing of the operator in atomic contributions Finally, y 5 








is one of the Dirac matrices, Eq. ([3]), with 12x2 and 02x2 being 
the 2 x 2 identity and null matrices, respectively 

02x2 hx2 /- 

.12x2 0 2x2 J ' W 

The Fermi coupling constant Gp = 2.22254 x 10 _14 ^h«o 
demonstrates the minuteness of the PV effect. 

We have calculated PV energies at the HF level as well as 
using the density functionals B3LYF®^and PBE,® all in a 
4-component relativistic framework. For the DC single-point 
calculations we have employed uncontracted cc-pVDZ basis 
sets SI for all atoms except Re, for which we have used the 
basis sets employed in Ref. 25 The small component ba¬ 
sis sets were generated by restricted kinetic balance imposed 
in the canonical orthonormalization step .El The two-electron 
Coulomb integrals (SS|SS), involving only the small compo¬ 
nents, were neglected in all calculations and the energy cor¬ 
rected by a simple point-charge model . 55 
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